Phase Diagram and Calorimetric Properties of NaFei_a;Co x As 
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We measured the resistivity and magnetic susceptibility to map out the phase diagram of single 
crystalline NaFei-^Coa; As. Replacement of Fe by Co suppresses both the structural and magnetic 
transition, while enhances the superconducting transition temperature (T c ) and superconducting 
component fraction. Magnetic susceptibility exhibits temperature-linear dependence in the high 
temperatures up to 500 K for all the superconducting samples, but such behavior suddenly breaks 
down for the non-superconducting overdoped crystal, suggesting that the superconductivity is closely 
related to the T-linear dependence of susceptibility. Analysis on the superconducting-state specific 
heat for the optimally doped crystal provides strong evidence for a two-band s-wave order parameter 
with gap amplitudes of Ai(0)/fcBT c = 1.78 and A2(0)/£;bT'c=3.11, being consistent with the nodeless 
gap symmetry revealed by angle-resolved photoemission spectroscopy experiment. 

PACS numbers: 74.70.Xa,74.25.Dw,74.25.Bt,74.20.Rp 
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I. INTRODUCTION 

Superconductivity in the iron-based superconductors, 
just like that in cuprates, emerges in the proximity to 
magnetically ordered state, so that the magnetic inter- 
actions are considered to play a key role in the mech- 
anism of such high-T c superconductivity. Accordingly, 
one of central issues for the iron-based superconductors is 
whether SDW follows a local moment or a Fermi surface 
nesting picture because this is correlated to the nature of 
the pairing force responsible for the superconductivity*^— 
The interplay between magnetism and superconductiv- 
ity has been extensively investigated in doped AeFe2As2 
(Ae = Ca, Sr, Ba Eu, and K, so called "122") family due 
to easiness of obtaining sizeable and high-quality single 
crystals- 5 Phase diagram, through measurements of elec- 
trical transport, magnetism (susceptibility, fj.SK, Neu- 
tron scattering, NMR), crystal structure and so on, has 
been well studied in such 122 single crystals^— While 
for the doped NaFeAs system (so called "111"), stud- 
ies for phase diagram were performed only on polycrys- 
talline samples but not on single crystalline ones due to 
the hardness to growing high-quality single crystals and 
the difficulty to control the doping concentration. Ac- 
cording to the same reason, little work has been done 
to elucidate the symmetry of the superconducting gap of 
doped NaFeAs system ji^rilL although quite a few studies 
have been carried out to explore the superconducting gap 
structure in the isostructural LiFeAs by either specific 
heat or heat transport or ARPES. 16-20 Even for these 
few studies on the symmetry of the supercondcuting gap 
of doped NaFeAs, the results are inconsistent with each 
other. ARPES and heat transport experiments revealed 
nodeless gap s 14 ' 15 while the measurements of penetration 
depth suggested gaps with nodes^ 

NaFeAs is established in Fe^As structure with the in- 



terstitial Fe replaced by Na atoms. Therefore, NaFeAs 
consists of a common building block, the FeAs layer, and 
the double layer of Na + sandwiched between the FeAs 
layers. NaFeAs is expected to be a simplified version of 
the structure of -ReFeAsO (Re = rare earth, so called 
"1111") and "122" iron pnictides. There is no static 
magnetic ordering and structural transition in isostruc- 
tural compound LiFeA s 21 i 22 , which shows superconduc- 
tivity at 18 K. However, NaFeAs itself has a spin-density- 
wave (SDW) magnetism with a small magnetic moment 
(0.09 /xe/Fe) 23 , in contrast to the larger values of mag- 
netic moment ~ 0.4 ^b/Fs in the La-" 1111"- and ~ 0.9 
^s/Fe in the Ba-"122" parent compounds 24 . Though 
NaFeAs shows superconductivity without purposely dop- 
ing, it possesses only 10% superconducting volume frac- 
tion and long range AFM order exists in most of its vol- 
ume inferred from the fiSK study on the polycrystalline 
sample^. By doping with Co on Fe site, bulk super- 
conductivity in NaFeAs can be achieved. Systematical 
study on the doped NaFeAs single crystals by varying 
Co concentration has not been done yet. In this pa- 
per, we report the study on the phase diagram of sin- 
gle crystalline NaFei_ x Co K As by measuring resistivity, 
magnetic susceptibility and specific heat. Additionally, 
susceptibility shows temperature-linear dependence up 
to 500 K for all the superconducting samples, and the 
deviation from high-temperature linear behavior occurs 
in the low temperatures, and the deviation temperature 
increases slightly with increasing the Co concentration. 
However, such behavior is suddenly changed for the non- 
superconducting overdoped crystal. Analysis on the spe- 
cific heat of the optimally doped sample revealed a two- 
band s-wave order symmetry with two gaps with their 
sizes at T = to be Ai(0)/fc B T c = 1.78 and A 2 (0)/fc B T c 
= 3.11, suggesting the strong-coupling superconductiv- 
ity. 



* Corresponding author; Electronic address: chenxh@ustc.edu.cn 
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II. EXPERIMENTAL DETAILS 



High-quality single crystals of NaFei-^Co^As have 
been grown by a NaAs flux method. NaAs was obtained 
by reacting the mixture of the elemental Na and As in 
evacuated quartz tubes at 200 °C for 10 h. Then Fe, 
Co and NaAs powder were carefully weighed according 
to the ratio of NaAs: Fe: Co=4: 1-x: x with (x = 
- 0.3), and thoroughly ground. The mixtures were put 
into an alumina crucible and then sealed in an iron cru- 
cible under 1.5 atm of highly pure argon gas. The sealed 
crucibles were heated to 950 °C at a rate of 60 °C/h 
in the tube furnace filled with the inert atmosphere and 
kept at 950°C for 10 h, and then cooled slowly to 600 
°C at 3°C/h to grow single crystals. The shiny crys- 
tals with typical size of 5x5x0.2 mm 3 can be easily 
cleaved from the melt. X-ray diffraction was perform on 
Smartlab-9 diffracmeter (Rikagu) from 10° to 60° with 
a scanning rate of 2° per minute. The actual chemical 
composition of the single crystals is determined by en- 
ergy dispersive X-ray spectroscopy (EDS). With chang- 
ing the nominal x from to 0.3, we obtained 12 batches 
of crystals, whose actual compositions were determined 
by EDS to be 0, 0.006, 0.010, 0.014, 0.017, 0.021, 0.028, 
0.042, 0.047, 0.061, 0.070, 0.075 and 0.109, with stan- 
dard instrument error as about 10%. Measurements of 
resistivity and specific heat were carried out by using 
the PPMS-9T (Quantum Design). Magnetic susceptibil- 
ity was measured by using SQUID-MPMS-7T (Quantum 
Design), and a high-temperature oven was used in the 
SQUID-MPMS for magnetic susceptibility measurement 
above 400 K. In all the magnetic measurements, magnetic 
field was applied within ab-plane. 
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FIG. 1: (a): Selected XRD patterns for NaFei-zCoz As single 
crystals, (b): The lattice parameter of c-axis plotted as a 
function of Co concentration x. Inset in (a) is the rocking 
curve for the crystal with x — 0.028. 
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FIG. 2: Resistivity as a function of temperature for 
NaFei-zCoa; As single crystals. Inset is the derivative resistiv- 
ity for the underdoped NaFei-xCo^As single crystals, where 
the scenarios to determine the T s and Tsdw are given. 



III. RESULTS AND DISCUSSION 

Figure 1(a) shows the selected single- crystalline XRD 
patterns for NaFei-^Co^As single crystals. Only (00Z) 
reflections can be recognized, and the rocking curve plot- 
ted in the inset of Fig. 1(a) shows a full- width- at-half- 
maximum (FWHM) about 0.05°, indicating c-axis ori- 
entation and high quality for these single crystals. Lat- 
tice parameter of c-axis (c) estimated from Fig. 1(a) 
was plotted as a function of Co concentration, as shown 
Fig. 1(b). The lattice parameter, c, of the undoped com- 
pound is 7.056 A, being consistent with previous reported 
results ] 25 ! 26 The lattice parameter decreases linearly with 
increasing Co doping, consistent with the results reported 
on polycrystalline NaFei-^Co^As samples^ While the 
amplitude of the change of c (~ 1% from to 0.109) is 
much smaller than that reported in polycrystalline sam- 
ples (~ 4%) ,2* 

Temperature dependence of resistivity taken from 2 to 
300 K were displayed in Fig. 2 for the selected single- 
crystalline NaFei-zCo^As samples. The resistivity de- 
creases with increasing Co doping level. An upturn in 
resistivity is observed in the low temperatures for the un- 
derdoped crystals, which arises from the structural/ SDW 
transition. The inset shows the derivative of resistiv- 
ity to figure out the temperatures corresponding to the 
structural and SDW transition (T s and Tsdw)- The two 
distinct features in dp(T)/dT are used to determine the 
T s and Tsdw- The derivative of resistivity is shown in 
the inset of Fig. 2. The structural and SDW transitions 
are suppressed rapidly with increasing Co concentration. 
In samples with Co concentration higher than 2.1%, no 



3 



trace of structural/ SDW transition can be recognized in 
resistivity. The undoped compound is already super- 
conducting although the superconducting transition is 
quite broad, and the resistivity of NaFeAs reaches zero 
at around 10 K, consistent with previous reports. The 
superconducting transitions for most of the samples are 
quite round, so we define the temperature at which re- 
sistivity reaches zero as T c . The optimal doping level 
is around x — 0.028 and the corresponding T c is 20 K. 
Further Co doping leads to the decrease of T c , and no 
superconducting transition can be observed down to 2 K 
for the crystal with x = 0.109. 
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FIG. 3: Magnetic susceptibility taken at 10 Oe in ZFC mode 
for the superconducting samples. 

Figure 3 shows the magnetic susceptibility taken at 
10 Oe in zero field cooling (ZFC) procedure for the su- 
perconducting NaFe^a-COzAs single crystals. For x — 
0, tiny diamagnetic signal was observed below 9 K, in- 
dicating the superconductive shielding effect is weak in 
this composition of crystal. With enhancing Co doping 
level, T c inferred from diamagnetic signal increases. The 
T c determined by susceptibility is nearly the same as T c 
obtained from resistivity measurements for the crystals 
with the same composition. The superconductive shield- 
ing fraction rises steeply with increasing Co doping level. 
Full shielding at 2 K can be observed for the crystal by 
substituting for Fe with 1% Co. The highest T c of 20 K 
is obtained in the crystal with x — 0.028, and the super- 
conducting transition is very steep. To our knowledge, 
this is the best for all the reported crystals in this sys- 
tem. Both shielding fraction and T c decrease with further 
doping of Co. 

Susceptibility up to 500 K taken under the applied 
magnetic field of 5 T is shown in Fig. 4(a). Because of 
the very close magnitude of the normal-state susceptibili- 
ties for the superconducting samples, the susceptibilities 
are normalized to the values at 300 K and a shift was 
then made for all the Co doped samples to distinguish 
from each other, as shown in Fig. 4(b). Rapid drops due 
to superconducting transition can still be observed at low 
temperature for the doped superconducting samples. For 




FIG. 4: (a): Normal-state susceptibility measured at 5 T in 
the temperature range up to 500 K for the NaFei-zCoa; As 
single crystals; (b): Normalized normal-state susceptibility 
for the NaFei-^Coz, As single crystals. The susceptibility was 
normalized to the value at 300 K and shift from sample to 
sample to distinguish from each other. 



the undoped and underdoped samples, slight kinks can be 
observed in magnetic susceptibility just above T c between 
20 K and 55 K, which arises from the structural and SDW 
transitions. Temperature linear dependence can be ob- 
served in high temperatures for the magnetic susceptibil- 
ity of all the superconducting samples. Magnetic suscep- 
tibility deviates from the high-temperature linear behav- 
ior in the low temperatures, and the deviation tempera- 
ture lies at about 200 K for all the superconducting crys- 
tals. While with further increasing the Co doping level, 
only very weak linear behavior can be found above 400 
K in heavily overdoped and non-superconducting crystal 
with x = 0.109. The slope for the linear dependence of 
high-temperature susceptibility is the nearly the same for 
all the superconducting crystals. Such temperature lin- 
ear dependence of magnetic susceptibility is an universal 
feature for all the iron-based superconductors, which has 
been theoretically ascribed to the spin fluctuations aris- 
ing from the local SDW correlation^ Another theoreti- 
cal interpretation for the temperature linear dependence 
of the susceptibility is taking into account the peculiar- 
ities of orbitally resolved densities of states due to local 
correlations, which does not need to invoke the antifcrro- 
magnetic fluctuations =§, Whether the temperature linear 
dependent susceptibility in high temperatures is corre- 
lated with the pairing force in the high-T c superconduc- 
tivity of the iron-based superconducting compounds is 
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still open question. In the present system, breakdown 
of the linear dependent susceptibility in the overdoped 
region coinciding with disappearance of superconductiv- 
ity seems to suggest the association between temperature 
linear dependence of the susceptibility and superconduc- 
tivity. 
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FIG. 5: Phase diagram of NaFei-xCoi As. T s , Tsdw and T c 
are determined from the resistivity in Fig. 2 and its inset, 
susceptibility in Fig. 3 and Fig. 4 give almost the same value 
for T s , Tsdw and T c . 

Based on Ts, Tsdw an d T c inferred from the resis- 
tivity and magnetic susceptibility, the phase diagram of 
NaFei-zCozAs system is established from the measure- 
ments on single crystals, as shown in Fig. 5. The grad- 
ual destruction of magnetism and enhancement of su- 
perconductivity (rise of T c and superconductive shielding 
fraction) were observed with increasing Co doping level. 
A dome-shaped T c vs. x relationship can be observed. 
Optimal T c was obtained around x = 0.028 and further 
Co doping suppressed superconductivity. T c goes to zero 
around x = 0.109. This phase diagram is quite similar to 
those of 122 and 1111 systems except for that the starting 
compound of the present system is superconducting. In 
addition, the optimal T c can be achieved by about 2.8% 
Co doping in NaFei-^Co^As, much less than ^7% Co in 
Ba(Fei_ x Co x )2As2 system.] 29 ' 30 From this point of view, 
the phase diagram of NaFei_ x Co x As is more similar to 
that of Ba(Fei_ 2; Ni a ;)2As2 system^ in which the opti- 
mal doping is reached at 4.6%, far from that of ~7% in 
Ba(Fei^ x Co x ) 2 As2. 

Figure 6 displays the temperature dependence of the 
low-temperature specific heat (plotted as C p /T) for the 
underdoped and heavily overdoped crystals. The specific 
heat for the undoped samples shown in Fig. 6(a) ex- 
hibit two anomalies corresponding to the structural and 
SDW transitions, respectively. The Ts and Tsdw are 
well consistent with those determined from resistivity in 
the inset of Fig. 2. No anomaly corresponding to the 
superconducting transition can be observed in the spe- 
cific heat for the undoped sample, which could be due to 
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FIG. 6: (a): Temperature dependence of low-temperature 
specific heat (C p /T) for underdoped and heavily overdoped 
NaFei-zCoz As crystals. The arrows points to the anomalies 
in the specific heat, (b): The inset shows the derivative of 
specific heat for the crystal with x — 0.014, where structural, 
SDW and superconducting transitions can be clearly recog- 
nized from the dips, (c): Specific heat as a function of temper- 
ature for the optimally-doped and overdoped superconduct- 
ing NaFei-xCoa; As crystals. The blue line is the fitting of the 
specific heat between 20 - 40 K by using C' p = Cen+Ciatuce- 
(d): Electronic specific heat C e /T (by subtracting the lat- 
tice contribution from C v /T) as function of temperature for 
the optimally doped crystal, where the dashed line represents 
the normal-state electronic contribution, j n = 10.1 mj mol -1 
K- 2 . 
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the low fraction of superconducting component. This is 
consistent with the tiny superconducting shielding frac- 
tion from susceptibility measurement shown in Fig. 3. 
For the crystal with x — 0.014, a clear anomaly corre- 
sponding to the supercondcuting transition around 14 
K, indicating bulk superconductivity in this sample and 
consistent with the good superconductive shielding for 
crystals with x > 0.01 (as shown Fig. 3). In addition, 
an anomaly is observed at about 33 K, which coincides 
with the structural transition determined by the deriva- 
tive of resistivity shown in the inset of Fig. 2. However, 
no clear anomaly corresponding to SDW transition is ob- 
served. In order to distinguish from the structural, SDW 
and superconducting transitions, the derivative of spe- 
cific heat is shown in Fig. 6(b), and Fig. 6(b) clearly shows 
three dips at Tg, Tgnw and T c , indicating the coexistence 
of the antiferromagnetism and superconductivity for the 
underdoped crystal with x=0.014. The specific heat of 
the heavily overdoped and non-superconducting crystal 
with x = 0.109 does not show any anomaly. As temper- 
ature increases, the specific heat above T s for the sample 
with x = 0.014 tends to be the same as that of the sample 
with x = 0.109. 

Pronounced jump due to superconducting transition 
can be observed in the temperature dependence of spe- 
cific heat for the optimally-doped and overdoped super- 
conducting samples, as shown in Fig. 6(c). The normal- 



state specific heat, C p , can be described by C p 



+ Ciattice(T') with the electronic contribution of C en = 
7 n T and the lattice contribution of Ciatticcp 1 ) = (3T 3 
+ ST 5 . The solid line in Fig. 6(c) is the best fit to the 
C p /T above T c (20 K -40 K) for the optimally doped sam- 
ple, yielding 7„ = 10.1 mJ/mol K -2 , j3 — 0.23 mJ/mol 
K" 4 and S = -0.0589 //J/mol K~ 6 . After subtracting this 
Ciattice component, the obtained results are shown in Fig. 
6(d), for which the entropy conservation is confirmed to 
be satisfied. Figure 6(d) shows that C e /T does not ex- 
trapolate to zero at T = limit but to a finite residual 
normal-state-like contribution 7,. = 1.5 mJ mol -1 K~ 2 . 
Finite values of 7 r are a common feature in the specific 
measurements of iron-based superconductors^ 

Normalized electron specific heat is displayed in Fig. 
7, where -f es = C e /T. Figure 7(a) indicates that C e can- 
not be described by a single band BCS superconductor, 
calculated either in the weak-coupling limit with a — 
A(0)/fc B T c = 1.764 or by letting a adjustable. Espe- 
cially, the fitting in weak-coupling limit shows very poor 
agreement with the experimental data. Figure 7(a) also 
shows the fitting results by using the single-band d-wave 
model. It is clear that such a fc-dependent gap, even 
in strong coupling scenario (obtained a = 3.09), cannot 
describe the observed data. As a result, the phenomeno- 
logical two-band a model with two energy gaps was used, 
and the fitting results are shown in Fig. 7(b). The two- 
band a model results in an excellent fitting of the specific 
heat from T c /10 to T c , and as a consequence gives reli- 
able gaps. This fitting is calculated as the sum of the 
contributions from two bands by assuming independent 
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FIG. 7: (a): The normalized electron specific heat of the opti- 
mally doped crystal with £=0.028, compared with the specific 
heat of single-band s-wave (wine line) and d-wave (magenta 
line) order parameters as well as that in the weak-coupling 
limit (red line), (b): The normalized electron specific heat of 
the superconducting sample with £=0.028. The blue curve 
represents a two-gap fit. The wine and green curves are the 
partial specific-heat contributions of the two bands. 



BCS temperature dependences of the two superconduct- 
ing gaps, as shown by the dashed and short-dashed lines 
in Fig. 7(b). In the fitting, the magnitudes of two gaps 
at T — limit are introduced as adjustable parameters, 
ai = Ai(0)/fc B T c and a 2 = A 2 (0)/A:bT c . At the same 
time, 7i/7n (* = 1, 2), which measures the fraction of 
the total normal electron density of states, are intro- 
duced as adjustable parameters. The obtained a\ for 
small gap is 1.78, and «2 for large gap is 3.11. The rel- 
ative weight for the small and large gaps here is 71/72 
= 0.72. Both of the obtained two gaps are larger than 
that of the BCS weak-coupling limit. This is not consis- 
tent with the theoretical constraints that one gap must 
be larger than the BCS gap and one smaller in a weakly 
coupled two-band superconductor, 32 which has been ob- 
served in MgB2^ and another iron-based superconductor 
Ba(Feo.925Coo.o75)2As2 — and as well as another 111 su- 
perconductor LiFeAs.— Actually, such phenomenon that 
the two gaps are larger than the BCS weak-coupling 
limit has been also observed in Fe(Teo. 57860.43) — and 
Bao.6Ko.4Fe2As2 — single crystals. The large mag- 
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nitudes of gaps derived from our fitting suggest the 
strong-coupling superconductivity, as that reported in 
Fe(Teo.57Seo.43) previously.— The relative ratio of the 
two gaps in NaFeo.972Coo.028 is Ai(0)/A 2 (0) ~ 0.57, 
which is a little larger than that seen in other iron- 
pnictide superconductors where Ai(0)/A2(0) = 0.3 ~ 
q 5 1 34 1 37- - 39 j^. should De pointed out that although in the 
present work the fitting of specific heat suggests isotropic 
gaps, some theoretical and other experiments also suggest 
complicated pair symmetry for most of iron-based super- 
conductors, including the nodeless anisotropic or nodal 
gaps. Actually, penetration depth experiment indicates a 
nodal gap in the optimally doped crystal^ while ARPES 
and heat transport measurements revealed nodeless gaps 
in the optimally doped and overdoped NaFei_ x Co x As 
crystal ) ' 15 Our results supports the nodeless gap sym- 
metry observed by ARPES experiment.— Especially, our 
analysis is consistent with the existence of different mag- 
nitudes of the gaps suggested by heat transport.— 

IV. CONCLUSION 

In summary, we mapped out the phase diagram of 
NaFei-zCozAs system through measuring transport and 
magnetic properties on high-quality single crystals. Sub- 



stitution of Co for Fe destroys both the structural and 
magnetic transition, while enhances the superconducting 
transition temperature (T c ) and superconducting shield- 
ing fraction. Susceptibility exhibits temperature-linear 
dependence in the high temperatures up to 500 K for all 
the superconducting samples, while such behavior cannot 
be observed in the non-superconducting overdoped crys- 
tal. The breakdown of the linear dependent susceptibility 
in the overdoped region coinciding with disappearance of 
superconductivity seems to suggest the association be- 
tween temperature linear dependence of the susceptibil- 
ity and superconductivity. Analysis on the specific heat 
in superconducting state reveals a two-band s-wave order 
parameter with gap amplitudes of Ai(0)/fcsT c = 1.78 and 
A 2 (0)//sbT c =3.II, suggesting the strong-coupling super- 
conductivity. 
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